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Abstract 

Applying a technique developed recently [1,2] for an harmonic oscillator coupled to a 
bath of harmonic oscillators, we present an exact solution for the tunneling problem in an 
Ohmic dissipative system with inverted harmonic potential. The result shows that while 
the dissipation tends to suppress the tunneling, the Brownian motion tends to enhance the 
tunneling. Whether the tunneling rate increases or not would then depend on the initial 
conditions. We give a specific formula to calculate the tunneling probability determined by 
various parameters and the initial conditions. 
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1. Introduction. 



Quantum tunneling with dissipation has been studied by many people after the work 
by Caldeira and Leggett [3-6] . These studies use different approximations to calculate the 
tunneling probability. Among these works, there is a widely discussed question about whether 
the dissipation suppresses or enhances the quantum tunneling[7]. In this paper, we shall show 
that for a special potential barrier, the inverted harmonic potential well, the tunneling problem 
can be solved exactly, thus we shall answer this question precisely. The technique we used 
for this solution has been developed recently to solve the wavefunction evolution for another 
dissipative system problem: an harmonic oscillator coupled to an environment [1,2]. In this 
paper we shall show that this technique is equally applied to the inverted harmonic potential. 

The simplest example of a dissipative system, an harmonic oscillator coupled to an envi- 
ronment of a bath of harmonic oscillators, has been the subject of extensive studies (see [1] 
for references). In a recent paper [1], we obtained a simple and exact solution for the wave 
function of the system plus the bath, in the Ohmic case. It is described by the direct product 
in two independent Hilbert spaces. One of them is described by the Caldirola-Kanai (CK) 
Hamiltonian, the other represents the effect of the bath, i.e., the Brownian motion. In a sec- 
ond paper [2], using this wavefunction solution, we derived a simple formula for the probability 
distribution at finite temperature, expressed in terms of the wavefunction solution for the CK 
Hamiltonian only. 

In this paper we shall apply these results to a different system, the inverted harmonic 
potential well to study the quantum tunneling problem of dissipative system. We shall use 
the same notations as the reference [1,2] and quote formulas from there. The Hamiltonian of 
this system is: 
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The only difference of this Hamiltonian from that in reference [1] is the sign of the second 
term of the right hand side, i.e., the potential is now an inverted harmonic potential well. 
Assuming the Ohmic bath oscillator density to be the same as given in [1]: 

= —-f- (2) 

Following exactly the same procedure as [1], we derive the Langevin equation of motion for 
the system: 

q + rjq - ufa = /(f) , (3) 
with the Brownian motion driving force: 

r / \ X — ^ ^ 7 / • ^ ^ 7 ^ \ / a \ 

/(*) = - E -6(x j0 cos ujjt + x j0 -^-) . (4) 

This equation is easily solved to give: 

q(t) = ai(t)g + a2<»<?o + E (MO^jo + b j2 (t)x j0 ) , (5) 

j 

where <?o, 3o, ^jO) ^jo are the initial position and velocity operators of the system and bath 
respectively, and 

-it, n V ■ i \ _v.fSmh.ujt 
ai — e 2 (coshcjt + — smhut), a 2 = e ? , (6) 

b j± (t) = -|l J* a 2 (0 coswj (f - df', 6j-i (f) = -jj J* a 2 (t') sinujj (t - t') dt', (7) 

with u = (ui + ^y/ 2 . 

Then, based on the same arguments of [1], this solution can be used to show that the 
wavefunction of the system plus the bath can be written as 

*(?,{&}>*) = #z-Efc>*) lixA^t) , (8) 



where the wavefunction if)(Q, t) is a solution of Schoedinger equation with the time dependent 
CK Hamiltonian 

H Q = eT nt Mule^Q 2 , (9) 

and the commutation [Q,P] = ih, while ^ = bji(t)xj + bj2(t)xj is the contribution of the 
j'th bath oscillator to the Brownian motion. The function Xj{£jit) i s given by Xji^j^t) = 
< 9^ \xjo >, where Xjo( x jo) is the initial state of the bath oscillator, and \9^. > is an eigenstate 
of the operator £ 3 - with eigenvalue as given by eq.(4) of [2]: 



V'; l-''/U ~ -''/t's./ 



'j2 



(10) 



According the analysis of [2], assuming the bath is at temperature T, averaging over the 
Boltzmann distribution of the wavefunction Xjo{ x jo) of the bath oscillators, and using eq.(§), 
we find the probability distribution the same as derived in [2]: 



c2 



pUi.t) = -^~ I \rUi-^.t)\ 2 e 2 V c /(. (11) 



where is the Brownian motion width: 



crf(t) = 1 l^ihM' + ^W)?) coth ^) piu^. (12) 
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This result is apparently the same as that given in [2], but the expressions for bj\(t) and 
bj2 (t) are different. There are also subtleties associated with a logarithmic divergence of the 
integration over Uj, which is removed by introducing a cut-off frequency u c . This width is 
zero initially, but then increases exponentially to infinity as time evolves. In the following, for 
simplicity, we consider only the low temperature limit, then the width approaches in a time 
range of about 1/u to the asymptotic value 

4 = e^K^B, (13) 



with g 2 = h/ (2Mujq) (the width of the ground state of an harmonic oscillator with frequency 
ojq) and 
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B^^W + i-^n (14) 
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2. The wavefunction evolution 

To calculate the tunneling probability, we assume that initially the wavefunction is a 
gaussian wave packet centered at the right of the peak of the potential by z with a velocity 
v = —hk/M towards the potential barrier peak: 

(g ~ z ^ | ika 

MQ) = t = 0) = (2na 2 y 1 / A e 4o* \ (15) 

First, following exactly the same method in [1], we derive the Green's function of the time 
dependent Hamiltonian eq.@: 



, / M \ 2 r iM , 2 _* 2 
G(g,g ;^,0) = y— —J exp [a^ + a 2 e{ q - 2q q 

Then we calculate the wavefunction ip(q,t) : 



(16) 



ip(q,t) = J G(q,q ;t,0)ip (q )dq (17) 



(q ~ z ) 2 



= (27ra 2 )- 1/4 (a 1 + iu a 2 r 2 ) e 4<x 2 (18) 

This result represents a gaussian distribution centered at q c = aiz — a 2 v ( the classical 
trajectory of a particle initially at zq with velocity -Vq) and with a width of 

a 2 = a 2 (a 2 +r 4^2 a 2 )5 (lg) 

where r = gq/g. The coefficients in the phase factor exponent are: 



C2 = IFa 1 """'' (20) 

ci = -^-(-2g c -(lna 2 ) + 4g c ), (21) 

co = ^(g c ^(lna 2 ) - 2g c ) + f§o; Q a 2 r 2 . (22) 
4 at 4o"g 

This phase factor is related to the current density, as will be explained later. It is straight 
forward to verify that eq.(|18D indeed satisfies the Schoedinger equation with the time de- 
pendent Hamiltonian eq.(^). With these provisions, we are ready to calculate the tunneling 
probability. 

3. The tunneling probability and current density 



Using the density eq.(^TJ) and wavefunction eq.([Tj|), we find a very simple expression for 
the probability density: 

1 _(q-qc) 2 

p(q,t) = -^e ^T, (23) 

where 



o* = al + o% (24) 

is the total width including the Brownian motion width. 

The probability for the particle to pass to the left of position q at t is then: 



P(q,t)= I" p(q',t)dt' = F( q -^), (25) 



oo 



where the function F is: 



1 r w 2 

F(W) = -j= / e~ u du. (26) 
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Hence the current density is found to be: 



dP(q,t) |/12 .1. .d n 2 



-(g-g c )-(lna t 2 ) + g c ). (27) 



<9t Iri v 2^ ^dt 
On the other hand, if the Brownian motion can be ignored (replacing at by oq ) , the solution 
eq.flT8|) of the time dependent Schoedinger equation should also satisfy the expression for the 
current density: 



The expressions for ci and c 2 eq.(|21"D and eq. (|20|) indeed satisfy eq.(p7|), hence they acquire 
the following clear physical meaning. If we denote the velocity of the "fluid" distribution |^| 2 
as v, the current density is then / = \ip\ 2 v. Therefore a comparison of eq. (|27D with eq.(^8|) 
leads to: 



h -(2c 2 q + c 1 ) = hq-q c )^-(\na 2 e ) + q c . (29) 



Thus those particles with larger q acquire larger velocity falling from the apex of the potential. 
Since the width a 2 e is proportional to exp(2c<j — r])t, as time evolves, most particles move rapidly 
away from the apex of the potential. In the mean time, for any fixed q ^ 0, the wavenumber 
2c 2 q + C\ also increases exponentially by a factor e vt according to eq. fl2~l~|) and eq. (pDD because 
the equivalent mass is Me vt instead of M. Hence the wavelength decays exponentially, as 
pointed out by reference [2]. When t 3> 1/rj the system approaches the classical limit. 

4. Tunneling probability as time approaches infinity 

We define the tunneling probability P as the probability for the particle to be at the left 
of the peak. Then, applying eq. fl25D , eq. ([24]) and eq.(|i~9~D, we find P = F(— W),mth 



V2a t ^2(a 2 (a 2 + r^u 2 a 2 ) + (r 2 ) 
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Because a±, 02, and are all proportional to e^'^l 2 ^ 1 as time approaches infinity, W 
approaches a finite limit, which determines the final tunneling probability. To simplify the 
expression for this limit, we use the following scaled variables: 

Z = — ; V= ; e = — ; r = — ; B = -In 1 + - . 31 

cr ujqZq 2uj a 71 U ~ 2 

Then, using the eq.(|6]) for a\a2 , and eq.(|T3|) for the Brownian motion width, we have: 

W = 4= , V (32) 

V2 ^1 + i^fr 4 + ABeVT^r 2 

To understand the meaning of this expression, we remark that Z ^> 1 represents a case 
where the wavepacket width is much smaller than its distance from the origin. We also remark 
that if 7] = 0, as time approaches infinity, q c = e UJ ° t {zQ — vq/ujq)/2 . Hence V < 1 means that a 
classical particle with initial velocity vq and position zq does not have enough kinetic energy 
to pass the potential barrier if there is no dissipation. Thus Z ^> 1 and V < 1 represent a 
case relevant to the quantum tunneling problem. 

The term ABe^/T^ e 2 r 2 comes from the Brownian motion. If e 1, i.e., if the dissipation 
is small, this term increases as e increases, and reduces W, which in turn increases F(—W), 
hence enhances the tunneling. This effect is reduced if r <C 1, or, in other words, if the 
initial wavepacket width is much larger than o"othe increase due to The Brownian motion is 
insignificant. If the initial velocity is not zero, the second term in the numerator increase W 
as e increases, thus the damping suppresses the tunneling. Intuitively this can be explained 
as that the damping makes the classical particle unable to move to the barrier peak as close 
as if there is no damping, and hence increases the barrier height. 

To get an idea about the effect of dissipation on the tunneling probability, in figure 1, we 
plot P as a function of V and e for Z = 3, r = 0.3, and B = 3, which corresponds to a cut-off 
frequency about 100 time larger than u. Because the logarithmic dependence of B on u c , the 
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result is very insensitive to the cut-off frequency. The plot clearly shows that when V is large, 
the tunneling is suppressed by the dissipation, while if initial velocity is zero, the dissipation 
enhances the tunneling. 

5. Tunneling probability and uncertainty principle 

Finally, it is interesting to examine the relation between the initial momentum and position 
distribution and the tunneling probability. For simplicity, we consider here only the case 
without dissipation. The Fourier transform of initial wavefunction eq . (|15|) , the wavefunction 
in momentum space is: 

^J k ) = ^y/4 e -(k +k0 )i-i(k + k )z _ (33) 

7T 

Thus from the classical point of view, the probability for the particle to pass the potential 
barrier is the probability of the initial velocity v < —u z , or k < —z /2<Jq (notice that 
v = hk/M): 

P = jT^ \Mk)\ 2 dk = F(-^01). (34) 

A comparison with eq.([32|) shows that if r 3> 1 , i.e., if the initial width a is much smaller 
than the minimum wavepacket width o"o, or, in other words, if the initial momentum spread 
is very large, this crude estimate is correct. 

Another extreme is when r<l, i.e., when we can consider the initial momentum spread 
is very small, and all the particles have velocity t> , but initial position has large spread. Then 
from the classical point of view, the probability for particle to pass the potential barrier is its 
initial q < vq/ujq, i.e.: 



\Ml)\ 2 dq = F(- 1 V} ). (35) 

-oo V I 



Again, a comparison with eq.fl32|) shows that if r <C 1 , i.e., if the initial width a is much 
larger than the minimum wavepacket width a , this crude estimate is also correct. These esti- 
mates indicate that the tunneling probability has a very simple relation with the uncertainty 
principle, and we have a very simple way to estimate the tunneling probability. 
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Figure Captions 



Figure 1. The tunneling probability as a function of the scaled dissipation coefficient e 
and the scaled initial velocity V, assuming the scaled initial position Z = 3, the scaled inverse 
wavepacket size r = 0.3, and the scaled cut-off frequency B = 3. 
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